Hole-burning Diffusion Measurements in High Magnetic Field Gradients 
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We describe methods for the measurement of translational 
diffusion in very large static magnetic field gradients by NMR. 
The techniques use a "hole-burning" sequence that, with the 
use of fringe field gradients of 42 T/m, can image diffusion 
along one dimension on a submicron scale. Two varieties of 
this method are demonstrated, including a particularly efficient 
mode called the "hole-comb," in which multiple diffusion times 
comprising an entire diffusive evolution can be measured within 
the span of a single detected slice. The advantages and disad- 
vantages of these methods are discussed, as well as their poten- 
tial for addressing non-Fickian diffusion, diffusion in restricted 
media, and spatially inhomogeneous diffusion. 
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1. INTRODUCTION 

Field-gradient NMR is a widespread research tool. The 
common ingredient to any of its applications is the la- 
belling of space by the Larmor precession frequency of a 
nuclear species in the presence of a static magnetic field 
gradient G: 



loq (z) = jH Q (z) =jG-r 



[1] 



Here 7 is the gyromagnetic ration of the nucleus, ujo is 
the Larmor frequency, and r is the position within the 
sample. Two applications of this spatial dependence are 
(1) measurement of structure (imaging), and (2) mea- 
surement of motion (flow and diffusion). Diffusion mea- 
surements have been performed since the earliest days of 
NMR research (g, |, gj, and are among the most widely 
used image contrast effects in magnetic resonance imaging 
(MRJQ.gj 

One particular area of interest has been performing 
NMR experiments in large field gradients. Recent work on 
3-dimensional NMR microimaging with applied gradients 
of 50 T/m have successfully achieved voxel volume reso- 
lution of 40 femtoliters for biological cell imaging, (g, gj 
Another source of large gradients are those in the fringe 
fields of NMR magnets (g, gj. For common supercon- 
ducting magnets, such gradients are often in the range 
of 50 T/m; for high- field resistive magnets (Ho ~ 30 T) 



such as at the National High Magnetic Field Laboratory, 
fringe field gradients greater than 200 T/m exist. NMR 
superconducting magnet facilities with Maxwell pair de- 
sign exist for the purpose of generating gradients of or- 
der 200 T/m(g). Other gradient sources are those out- 
side of devices designed to probe exterior material, such as 
the NMR-MOUSE([7^), NMR well-logging tools Q 
or the magnetic resonance force microscope filS]). Alterna- 
tively, the internal magnetic field gradients within materi- 
als with inhomogeneous magnetic susceptibility have been 
used to study porous structure. |7^) A common issue 
to many of these cases is that the RF excitation pulses 
are "soft" and do not uniformly excite the sample, either 
due spatial inhomogeneity in the Hi field, the limited fre- 
quency bandwidth of a finite-duration pulse, or both. Re- 
cent studies have been performed to fully characterize the 
spin evolution in CPMG sequences of many "soft" pulses 
to correctly extract diffusive information ^ldj, Other 
analyses have been performed to adapt multiple-quantum 
coherence sequences to the field gradient regime (jl^j. New 
techniques using "nutation echoes" formed through a com- 
bination of inhomogeneities in the static (Ho) and RF(i?i) 
fields have been developed (fli|, \2(\), and were included in 
a scheme showing the successful recovery of full chemical 
shift information in the presence of a static field gradi- 
ent of 50 mT/m.(fjl|, [j^j As these studies have shown, 
prospects and applications for field-gradient NMR capa- 
bility are growing. Thus, it is essential to adapt existing 
NMR techniques to the field gradient regime, as well as 
recognize capabilities that only large gradients provide. 

In this article, we describe methods for the measure- 
ment of translational diffusion in large static field gradi- 
ents in the fringe field of NMR magnets. These methods 
are of the "hole-burning" variety, in which long, low power 
RF pulses are used for spectrally (and thus spatially) se- 
lective irradiations prior to detection. The time evolution 
of such "holes" can be analyzed to extract diffusion infor- 
mation on a sub-micron scale. Such a selective excitation 
technique has been successfully applied in the past using 
internal magnetic field gradients to study porous structure 
with liquids and using applied gradients to image the 
diffusion of gases (25); the constrast of the present study 
with that work is in the geometry and spatial resolution 
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of the resonant slice with a much larger magnitude of the 
applied gradient. In addition to providing a viable alterna- 
tive to dephasing methods, the methods we describe have 
potential for deeper study of non-standard diffusion pro- 
cesses, such as anisotropic or restricted diffusion. 



2. HOLE-BURNING DIFFUSION SEQUENCES 

Two hole-burning sequences are sketched in Figj^. 
The left panel shows a standard hole-burning sequence, 
employed in many NMR experiments as well as other 
spectroscopies. [g^, ]2ty . Such an experiment con- 

sists of applying a long, low power pulse to irradiate a 
narrow band in frequency, given by Av « l/t p . In the 

presence of a static magnetic field gradient G, such a pulse 
irradiates a narrow spatial hole perpendicular to the gra- 
dient direction. For example, in a gradient of G=42 T/m, 
the thickness of a 1 H hole irradiated by a t p =l ms pulse is 
approximately Az = ^ = w 0.6 /im. After a diffu- 



sion period r has elapsed, a broad detection is performed of 
a large slice whose thickness is typically a few hundred /xm. 
In this case, this takes the form of a Hahn echo sequence 
with "hard" RF pulses. If r << T\, the spins irradiated 
with the burn pulse are "edited out," and do not appear in 
the broadly detected signal. As a function of the evolution 
time t, the labelled spins diffuse, widening the hole shape 
while conserving its area (in the absence of relaxation). 
This time evolution can be analyzed to extract a diffu- 
sion coefficient. This technique works in practice because 
with large gradients the hole thickness can be brought to 
the scale of the diffusion length for an NMR experiment, 
which is typically a few /im for liquids. 

In the single-hole sequence, we must wait several spin- 
lattice relaxation times (Xi) after each acquisition for the 
magnetization to return to equilibrium. For long T%, ac- 
quiring spectra at many values of the evolution time t 
is time-consuming. A variation on the hole-burning se- 
quence, sketched on the right in Fig. [I], circumvents this in- 
convenience by using more of the available detection slice. 
In this sequence, not one but a series of hole-burn excita- 
tion pulses are applied. They are spaced out in time on 
the scale of the diffusion time, and each is applied at a 
different frequency (i.e. position) within the bandwidth of 
the broad detection pulse. After this "hole-comb" , the en- 
tire slice is detected, with the result sketched in the lower 
right panel of Fig. |l|. The spins in the earlier holes diffuse 
while later holes are burned, so that the final spectrum 
contains a set of snapshots comprising an entire hole evo- 
lution. This technique assumes a uniform diffusion coeffi- 
cient and field gradient across the slice, so all holes broaden 
at the same rate; for bulk liquids in fringe-field gradients, 
this uniformity is excellent (within 0.1%). The evolution is 
acquired in a single transient; consequently, this sequence 



has a dramatically improved efficiency, and therefore pro- 
vides higher sensitivity. The savings in acquisition time 
provides a signal-to-noise enhancement of y/N, where N is 
the number of holes (or evolution times r) in the sequence. 
We note that the hole-comb measurements we performed 
were made possible through the fast frequency switching 
capability of the MagRes2000® spectrometer designed by 
A. P. Reyes. (f|) 

3. ANALYSIS 

In this section we describe briefly the procedures used 
to extract a diffusion coefficient from a hole-burning exper- 
iment in a fixed field gradient. We denote the hole-burned 
absorption spectra at each evolution time t as F (lu,t) and 
the unburned spectrum as Fq(u>). The first step is to define 
the hole shape from the larger slice shape: 



A(u,t) = 1- 



F {u) 



[2] 



The hole shape A{uj,t) measures the profile of labelled 
spins along the gradient direction as a function of time. 
The hole's time evolution with time will be governed by the 
diffusive propagator P (u>, r; a/, 0). This gives the proba- 
bility that a given spin at frequency u>' at time t=0 will 
be found at frequency u at time t — r. For free, isotropic, 
diffusion, this 1-D propagator is well known: 
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[3] 



Given the propagator and the initial hole profile A (u', 0), 
the profile at a later time is found simply by the convolu- 
tion product, 



A (w, t) = fP {uj, r; cj', 0) A (u/, 0) dw' . 



[4] 



For example, we can calculate the explicit form of 
A (w, r) for normal diffusion and a gaussian hole shape. 
Given a gaussian as an initial hole shape, i.e. 



A g (tu,0) = e 3 



Eq. 



then becomes 



Ag (W, T) 



- (^) -«£ , 



[5] 



[G] 



\J\-kDt 

This gaussian integral can be easily performed by complet- 
ing the square and using J e~ ax dx — yj. The result is 



A g (uj,t) 



exp 



° 2 (r) 



[7] 
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where 

a 2 (r) = a 2 + A 1 2 G 2 Dt [8] 

We see that this result conserves the area of the hole, as 
expected; the amplitude of the gaussian decays just as the 
width is enlarged. Once the initial width has been mea- 
sured and fixed (through a spectrum measurement immedi- 
ately following the burn pulse), all subsequent spectra can 
each be fit with a single adjustable parameter a. The re- 
sulting square-widths can then be fit to linear time depen- 
dence to extract the diffusivity D. Thus far, in our exper- 
iments, we have used square-wave pulses which would be 
expected to burn sinc-function rather than gaussian holes. 
However, the gaussian is a convenient phcnomcnological 
form that, as shown in the experiments below, accurately 
describes the broadening of the hole. 

The resolution limits of this mode of diffusion mea- 
surement are determined by the spin relaxation times, T\ 
and T2. Spin-lattice relaxation causes the tag placed on 
the burned spins to evaporate, and fills in the hole with- 
out broadening. In our measurements, we separately mea- 
sure the spin-lattice relaxation time T\ and constrain the 
hole area to decay as e~ T / Tl . For sufficiently fast spin- 
lattice relaxation compared to the diffusion time, the hole 
broadening due to diffusion is indetectable. The minimum 
hole thickness is another bound on the experiment, deter- 
mined cither by the diffusivity D (fast limit), or by the 
spin-spin relaxation time T2 (slow limit). The hole thick- 
ness is controllable only if negligible diffusion or spin de- 
phasing takes place during the burn pulse length t p . The 
combined conditions place the following lower bounds on 
measurable diffusivities by this method, depending on the 
controlling factor in the minimum hole size (diffusion or 
spin-spin relaxation) (p?5[). 
For diffusion-limited holes. 



2 1 1 
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and for relaxation-limited holes 
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2 2 111 
■It,' GD> 2Y 1 T 2 



[10] 



Eq. [Io[] is similar to the limiting diffusivity accessible 
from stimulated echo methods ([|, An advantageous 

application of this method would be to a system with long 
T2, such as was found, with 1 H decoupling, in the natural 
abundance 13 C signal in glassy glycerol in a previous NMR 
hole-burning study in a homogeneous field (g). 



their gaussian fits are shown in Fig. [| The spin-lattice 
relaxation time T\ in this case is of the order of seconds, 
much longer than the evolution times of the experiment 
(up to 30 ms). The resulting diffusivity is D=(4.8 ± 0.1) 
10 _6 cm 2 /s, which compares well with a diffusion measure- 
ment by a more standard stimulated echo dephasing mea- 
surement with the same sample and gradient. 

Fig. |] shows a measurement by a 1 H NMR hole-comb 
sequence of the diffusivity of glycerol- 13 C2 at T=296 K. 
The corresponding individual hole shapes and their gaus- 
sian fits are shown in Fig. @. In these fits, the area of 
the hole was constrained to decay as e T i , with a spin- 
lattice relaxation time of Ti=112 ms measured separately 
by a standard saturation recovery sequence. The result- 
ing diffusivity is Z?=(2.86 ± 0.05) 10 _8 cm 2 /s, which again 
compares well with a stimulated echo diffusion measure- 
ment. 

Finally, we show in Fig. ^ a comparison of tempera- 
ture dependences of the glycerol- 13 C2 1 H diffusion coeffi- 
cient measured by two different methods: stimulated echo 
and hole-comb measurements. The agreement is good for 
higher temperatures, until the point that the burn pulse 
interval (sa 1 ms) approaches the spin-spin relaxation time, 
T2- At this point the hole-comb sequence becomes unreli- 
able. 

5. APPLICATIONS 

For normal diffusion, a convenient form can be derived 
for the hole evolution given any initial hole shape. Al- 
ternatively, if the diffusive propagator is non-Fickian, this 
measurement can serve to map out its behavior by convo- 
luting with a known excitation function. This capability 
distinguishes the hole-burning sequence from the dephas- 
ing sequences. Given a spectral evolution profile from a 
hole-burning sequence, higher order moments can be cal- 
culated which provide more information on the propaga- 
tor. Specifically, so long as the propagator is stationary, 
i.e. depends only on the difference of the observation times 



t = t — t' and positions Z = z 



it can be shown that 



any spatial moment of an evolved hole spectrum A(z,t), 
defined by 

M*\A(, -m - Jdzjz- (z)) n A(z,t) 
M * [A (Z ' T)] = JdzA(z,r) 

is simply the sum of the initial hole moment and that of 
the diffusive propagator: 



4. EXPERIMENT 

Fig. U shows a measurement by a 1 H NMR hole- 
burning sequence of the diffusivity of propylene carbonate 
at T=295 K. The corresponding individual hole shapes and 



M: [A (z, r)] = M: [A (z, 0)] + M n z [P (Z, r)] [12] 

Thus, to the extent that any function can be recon- 
structed with knowledge of its moments, this provides a 
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FIG. 2: Hole-burning diffusion measurement by 1 H NMR in propylene at T = 295 K in a gradient of G = 42 T/m. 
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FIG. 3: Fits of hole-burning diffusion measurement by 1 H NMR in propylene carbonate at T = 295 K in a gradient of G = 42 
T/m. 
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FIG. 4: Hole-comb diffusion measurement by 1 H NMR in glycerol at T = 296 K in a gradient of G — 42 T/m. 
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FIG. 6: X H NMR comparison of diffusivity measurements by stimulated echo and hole-comb techniques in the same sample of 
glycerol- 13 C2, with an applied gradient of G=42 T/m. The two techniques agree until the spin-spin relaxation time Ti is of 
order the burn pulse width hum; below this temperature translational diffusion cannot be correctly discerned by the hole-comb 
technique. The error bars are computed from the statistical accuracy of the fits, and are not shown if smaller than the symbol 



method for determining a non-Fickian propagator's spa- 
tial dependence. 

The hole-burning sequence is potentially useful in the 
study of porous media; if the initial hole size can be 
made less than the smallest confinement length scale, the 
crossover from free to restricted diffusion can be observed 
by steadily increasing the initial hole size. Such a variation 
is similar to that accomplished by variation of the first in- 
terval in a stimulated echo sequence. Finally, the spatial 
resolution of the hole-burning sequence is well-suited to 
problems of transport near surfaces; the diffusion coeffi- 
cient in a liquid can be inspected on a sub-micron scale at 
arbitrary distances from a solid-liquid or solid-gas inter- 
face. 



6. CONCLUSIONS 

We have demonstrated a hole-burning NMR diffusom- 
etry technique in large magnetic field gradients (G = 42 
T/m). The advantages this technique offers include spa- 
tial resolution and a higher descriptive capability for non- 
standard diffusion behavior. Future directions include ap- 
plication to porous media and diffusion near surfaces on 
the sub-micron scale. 
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